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ABSTRACT: We study the two-point Topological Charge Density Correlator (TCDC) in lattice 
QCD with two degenerate flavours of naive Wilson fermions and unimproved Wilson gauge action 
at two values of lattice spacings and different volumes, for a range of quark masses. Configurations 
are generated with DDHMC algorithm and smoothed with HYP smearing. In order to shed light 
on the mechanisms leading to the observed suppression of topological susceptibility with respect 
to the decreasing quark mass and decreasing volume, in this work, we carry out a detailed study 
of the two-point TCDC. We have shown that, (1) the TCDC is negative beyond a positive core 
and radius of the core shrinks as lattice spacing decreases, (2) as the volume decreases, the mag- 
nitude of the contact term and the radius of the positive core decrease and the magnitude of the 
negative peak increases resulting in the suppression of the topological susceptibility as the volume 
decreases, (3) the contact term and radius of the positive core decrease with decreasing quark mass 
at a given lattice spacing and the negative peak increases with decreasing quark mass resulting in 
the suppression of the topological susceptibility with decreasing quark mass, (4) increasing levels 
of smearing suppresses the contact term and the negative peak keeping the susceptibility intact and 
(5) both the contact term and the negative peak diverge in nonintegrable fashion as lattice spacing 
decreases. 
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1 Introduction 

As a consequence of the reflection positivity and the pseudoscalar nature of the relevant local 
operator in Euclidean quantum field theory, the two-point Topological Charge Density Correlator 
(TCDC) is negative at arbitrary non-zero distances [1]. In the continuum theory, close to the origin 
the two-point TCDC is negative and singular. From power counting, the singularity ~ — | _x; | 8 up to 
possible logarithms and hence is non-integrable. In order to obtain a positive and finite space-time 
integral (susceptibility), the TCDC should have a positive non-integrable singularity at the origin 



As the authors of Ref. [1] pointed out long time ago, this has non- trivial consequences for the 
derivation and interpretation of the Witten-Veneziano expression [3] for the r\' mass. The negativity 
of the TCDC also has non-trivial consequences related to the nature of topological charge structure 
in QCD vacuum [4]. 

The issues related to two-point TCDC are best studied in the theory rigorously formulated on 
a Euclidean lattice. However, the lattice theory defined by a particular action may not be reflection 
positive. Fortunately, this is not a concern for the Wilson fermion. The lattice operator for the 
topological charge density q(x) may extend over several lattice spacings, and thus for sufficiently 
small x, the continuum like behaviors are not expected. Nevertheless, continuum properties are 
expected to emerge as lattice spacings become smaller and smaller. Specifically, on a lattice with 
lattice spacing a, TCDC remains positive within a radius r c , which is expected to shrink to zero 
as a — > 0. The first investigation of lattice spacing dependence of the radius of the positive core 
and the negativity beyond the positive core of TCDC in lattice QCD was carried out in Ref. [5] in 
the context of overlap based topological charge density in quenched QCD. Later, similar study was 
carried out [6] for a variety of lattice QCD actions with and without quarks where discretization 
errors appear only at &{a 2 ). 

Flavour singlet axial Ward-Takahashi identity relates the topological susceptibility %, which 
is the four- volume integral of TCDC, to the chiral condensate in the chiral limit [7, 8]. As a 
consequence, % vanishes linearly in the quark mass in the chiral limit. Furthermore, at a given value 
of the quark mass, % is suppressed as volume decreases [8,9]. As part of an on-going program [10, 
11] to study the chiral properties of Wilson lattice QCD (unimproved fermion and gauge actions), 
recently, we have demonstrated the suppression of topological susceptibility with decreasing quark 
mass in the case of unimproved Wilson fermion and gauge action [12, 13] where, the suppression 
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Table 1. Lattice parameters and simulation statistics. Here block, N2, N tr j and T refers to HMC block, 
step number for the force F%, number of HMC trajectories and the Molecular Dynamics trajectory length 
respectively. 

of x with decreasing volume was also shown. In order to shed light on the mechanisms leading to 
these suppressions, in this work, we carry out a detailed study of the two-point TCDC. 

2 Measurements 

We have generated ensembles of gauge configurations by means of DDHMC algorithm [14] using 
unimproved Wilson fermion and gauge actions with n f = 2 mass degenerate quark flavours. At 
/3 = 5.6 the lattice volumes are 16 x 32, 24 3 x 48 and 32 3 x 64 and the renormalized quark mass 
ranges between 25 to 125 MeV (MS scheme at 2 GeV). At j8 = 5.8 the lattice volume is 32 3 x 64 
and the renormalized physical quark mass ranges from 15 to 75 MeV. The lattice spacings are 
determined using nucleon mass to pion mass ratio and Sommer method. These determinations 
agree for the value of Sommer parameter ro = 0.44 fm. The lattice spacings at jS = 5.6 and 5.8 are 
0.069 and 0.053 fm respectively. The number of thermalized configurations ranges from 7000 to 
14000 and the number of measured configurations ranges from 200 to 500. 
The topological susceptibility 

X = jd 4 xC(r) (2.1) 
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Figure 1. (left) Effect of smearing on C(r) at j3 = 5.8 and K = 0.15462 at lattice volume 32 3 x 64. (right) 
Effect of smearing on the % at j3 = 5.8 and K = 0.15462 at lattice volume 32 3 x 64 (taken from Ref. [13]). 

with the TCDC, 

C(r) = (q(x)q(0)}, r=\x\ (2.2) 

where q(x) is the topological charge density. For q(x), we use the lattice approximation developed 
for SU (2) by DeGrand, Hasenfratz and Kovacs [15], modified for SU (3) by Hasenfratz and Neiter 
[16] and implemented in the MILC code [17]. It uses ten link paths described by unit lattice vector 
displacements in the sequence {x,y,z, — y, — x,t,x, — t, — x, — z} and {x, y,z, — x,t, — z,x, — t, — x, — y} 
plus rotations and cyclic permutations. We used HYP smearing with optimized smearing coeffi- 
cients a = 0.75, 0t2 = 0.6 and a?, = 0.3 [18]. Unless otherwise stated we have used 3 smearing 
steps in all our calculations. 

3 Results 

In order to extract the topological charge density reliably on the lattice, using the algebraic def- 
inition, smearing of link field is essential. Smearing however smoothens out short distance sin- 
gularities. Excessive smearing may in fact wipe out the fine details of the singularity structure. 
Both the positive and negative contributions to % w& affected in this manner. This is illustrated in 
Fig. 1 (left) where we show the effect of 3, 5, 10 and 20 HYP smearing steps on C{r) at j8 = 5.8, 
K = 0.15462 and lattice volume 32 3 x 64. However the susceptibility is remarkably stable under 
smearing after three smearing steps as illustrated in Fig. 1 (right) (taken from Ref. [13]). 

As was already stated, from theoretical considerations, we expect suppression of % with de- 
creasing volume at a fixed quark mass. In Fig. 2 we present the finite volume dependence of the 
C(r) at j8 = 5.6 and fc = 0.158 at lattice volumes 16 3 x 32, 24 3 x 48 and 32 3 x 64. We find that as 
volume decreases, the magnitude of the contact term and radius of the positive core decrease and 
the magnitude of the negative peak increases resulting in the suppression of topological suscepti- 
bility as volume decreases. 
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Figure 2. Finite volume dependence of the C(r) at /3 = 5.6 and K" = 0.158 at lattice volumes 16 3 x 32, 
24 3 x 48 and 32 3 x 64. 
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Figure 3. Comparison of the radius of the positive core of C(r) at two different lattice spacings for compa- 
rable pion mass. Lattice volume is 32 3 x 64. 
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Figure 4. Extraction of r c (0)/a using the data at two lattice points shown in Fig. 3. 

In Fig. 3, we compare the radius of the positive core of C(r) at j8 = 5.6 and 5.8 for compa- 
rable pion masses in physical units. The lattice volume is 32 3 x 64. The figure clearly exhibits 
the shrinking of the radius of the positive core of C(r) in physical units as one approaches the 
continuum. 

To verify that the radius of the positive core of C(r) indeed vanishes in the continuum limit, 
in Ref. [6] the data was fitted with the function r c {a)/a = C + Ba 2 since all the fermion actions 
considered in that work are 0(a) improved. The constant C was found to be close to 2. In Fig. 4, 
we present the data for r c (a)/a versus a. In our case, however, G[a) lattice artifacts are present 
and we fit the data to the formula r c /a = C + Aa. We obtain C « 2.7 and A m 17. Note that the 
lattice definition of q(x) used in Ref. [6] is different from what is used in the present work, and the 
value of C may crucially depend upon the definition of q(x). We need data at more lattice spacings 
in order to extract the coefficient of an 0(a 2 ) term potentially present in the data. 

In order to understand the detailed mechanism behind the suppression of topological suscepti- 
bility with decreasing quark mass, we need to investigate the quark mass dependence of the various 
features of the C(r). In Fig. 5 we present the quark mass dependence of C(r) with emphasis on the 
positive region at j8 = 5.6 and lattice volume 24 3 x 48. The magnitude of the contact term C(0) is 
seen to decrease with decreasing quark mass. In Fig. 6 we present the quark mass dependence of 
C(r) with emphasis on the crossover from positive to the negative region of C(f) and the negative 
peak region at j8 = 5.6 and lattice volume 24 3 x 48. The radius of the positive core and the mag- 
nitude of the negative peak of C(r) are seen to decrease and increase respectively with decreasing 
quark mass. The features presented in Figs. 5 and 6 result in the suppression of the topological 
susceptibility with decreasing quark mass. MILC collaboration [19] has made a similar observation 
regarding the dependence of the negative peak on quark mass. In Fig. 7 we present the correspond- 
ing topological susceptibilities (j3 =5.6 and lattice volume 24 3 x 48) as a function of the quark 
mass which clearly shows the suppression as quark mass decreases. This figure includes K = 0.158 
in addition to the fc's presented in Figs. 5 and 6. 
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Figure 5. The quark mass dependence of C(r) with emphasis on the positive region at j5 =5.6 and lattice 
volume 24 3 x 48. 
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Figure 6. The quark mass dependence of C(r) with emphasis on the crossover from positive to the negative 
region of C(r) and the negative peak region at j8 =5.6 and lattice volume 24 3 x 48. 

In Fig. 8 we show the lattice spacing dependence of the contact term and the negative peak of 
C(r) at comparable pion mass in physical units for j8 = 5.6 and 5.8 and lattice volume 32 3 x 64. For 
comparison, the corresponding quantities for pure gauge lattice theory at jS = 6.0983 (a = 0.078fm) 
and lattice volume 24 3 x 48 are also shown. Both the contact term and the negative peak of C(r) 
increase with decreasing lattice spacing, in accordance with the expectation from the continuum 
theory. 

From the definition of the topological susceptibility 

% = I d 4 xC(r) = [ 27l 2 (r 3 )drC{r) 
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Figure 7. Topological susceptibility at /3 =5.6 and lattice volume 24 x 48 as a function of the quark mass 
mass. 
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Figure 8. Lattice spacing dependence of the contact term and the negative peak of C(r) at comparable pion 
mass for j3 = 5.6 and 5.8 and lattice volume 32 3 x 64. For comparison, the corresponding quantities for pure 
gauge lattice theory at j3 = 6.0983 (a = 0.078fm) and lattice volume 24 3 x 48 are also shown. 



it is useful to define [19] a local susceptibility 

X (r)= f 2rc 2 (r'V C(r') 
Jo 



(3.1) 



in order to exhibit the lattice spacing dependence more clearly. In Fig. 9 we plot %(r) versus r 
at two lattice spacings at comparable quark masses. The contribution from the positive part of 
C(r) results in a peak at short distance. This is followed by a decrease due to the negative part of 
C(r). As lattice spacing decreases, the contribution from the positive part increases resulting in the 
increase of the peak of %{r). 
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Figure 9. The function ^(r), defined in Eq. (3.1) as a function of r at j3 = 5.6 and 5.8 at comparable pion 
masses. 
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Figure 10. Contributions to the susceptibility from positive and negative parts of C(r) at j3 =5.6 and 5.8 at 
comparable pion masses. 



Define the contributions to the susceptibility from the positive and negative parts of C(r) as 

Xp = 2n 2 {r' 3 )dr C{r) and (3.2) 
Jo 

2K 2 (r' 3 )dr' C(r') (3.3) 

respectively. According to the expectations from continuum theory, the negative singularity close 
to the origin and the positive singularity at the origin are both nonintegrable. Thus the contributions 
to % from positive and negative parts of C(r) are expected to diverge, nevertheless resulting in a 
finite x due to cancellation. In Fig. 10, we plot the contributions to the susceptibility from positive 
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and negative parts of C(r) at j8 = 5.6 and 5.8 at comparable pion masses. The data exhibited in 
Fig. 10 are in accordance with these expectations. 

It is known that the topological susceptibility decreases with decreasing quark mass and de- 
creasing volume. To understand the mechanisms leading to these suppressions, in this work, we 
carry out a detailed study of the two-point TCDC. We have shown that, with naive Wilson fermion 
and gauge action, (1) the two-point TCDC is negative beyond a positive core and radius of the core 
shrinks as lattice spacing decreases, (2) as volume decreases, the magnitude of the contact term and 
the radius of the positive core decrease and the magnitude of the negative peak increases resulting 
in the suppression of topological susceptibility as volume decreases, (3) the contact term and ra- 
dius of the positive core decrease with decreasing quark mass at a given lattice spacing and the 
negative peak increases with decreasing quark mass resulting in the suppression of the topological 
susceptibility with decreasing quark mass, (4) increasing levels of smearing suppresses the contact 
term and the negative peak keeping the susceptibility intact and (5) both the contact term and the 
negative peak diverge in nonintegrable fashion as lattice spacing decreases. 

Acknowledgements 

Numerical calculations are carried out on Cray XD1 and Cray XT5 systems supported by the 
10th and 11th Five Year Plan Projects of the Theory Division, SINP under the DAE, Govt, of 
India. We thank Richard Chang for the prompt maintainance of the systems and the help in data 
management. This work was in part based on the public lattice gauge theory codes of the MILC 
collaboration [17] and Martin Liischer [14]. 

References 

[1] E. Seiler and I. O. Stamatescu, Some remarks on the Witten-Veneziano formula for the Eta-prime 
mass, MPI-PAE/PTh 10/87, unpublished; E. Seiler, Phys. Lett. B 525, 355 (2002). 

[2] M. Aguado and E. Seiler, Phys. Rev. D 72, 094502 (2005) [hep-lat/0503015]. 

[3] E. Witten, Nucl. Phys. B 156, 269 (1979); G. Veneziano, Nucl. Phys. B 159, 213 (1979). 

[4] I. Horvath, S. J. Dong, T. Draper, F. X. Lee, K. F. Liu, N. Mathur, H. B. Thacker and J. B. Zhang, 
Phys. Rev. D 68, 114505 (2003) [hep-lat/0302009]. 

[5] I. Horvath, A. Alexandra, J. B. Zhang, Y. Chen, S. J. Dong, T. Draper, K. F. Liu, N. Mathur, 
S. Tamhankar and H. B. Thacker, Phys. Lett. B 617, 49 (2005) [hep-lat/0504005]. 

[6] F. Bruckmann, F. Gruber, N. Cundy, A. Schafer and T. Lippert, Phys. Lett. B 707, 278 (2012) 
[arXiv: 1107.0897 [hep-lat]]. 

[7] R. J. Crewther, Phys. Lett. B 70, 349 (1977). 

[8] H. Leutwyler and A. V. Smilga, Phys. Rev. D 46, 5607 (1992). 

[9] For a detailed discussion, see S. Durr, Nucl. Phys. B 611, 281 (2001) [hep-lat/0103011]. 

[10] A. K. De, A. Harindranath and S. Mondal, Phys. Lett. B 682, 150 (2009) [arXiv:0910.561 1 [hep-lat]]. 

[11] A. K. De, A. Harindranath and S. Mondal, JHEP 1107, 117 (2011) [arXiv:l 105.0762 [hep-lat]]. 

[12] A. Chowdhury, A. K. De, S. De Sarkar, A. Harindranath, S. Mondal, A. Sarkar and J. Maiti, Phys. 
Lett. B 707, 228 (2012) [arXiv:l 1 10.6013 [hep-lat]]. 



-9- 



[13] A. Chowdhury, A. K. De, S. De Sarkar, A. Harindranath, S. Mondal, A. Sarkar and J. Maiti, PoS 
LATTICE 2011, 099 (2011) [arXivrl 11 1.1812 [hep-lat]]. 

[14] M. Luscher, Comput. Phys. Commun. 156, 209-220 (2004). [hep-lat/03 10048]; M. Liischer, Comput. 
Phys. Commun. 165, 199-220 (2005). [hep-lat/0409106]. 
http : / /luscher . web . cern . ch/luscher/DD-HMC/ index . html 

[15] T. A. DeGrand, A. Hasenfratz, T. G. Kovacs, Nucl. Phys. B505, 417-441 (1997). 
[arXiv:hep-lat/9705009 [hep-lat]]. 

[16] A. Hasenfratz, C. Nieter, Phys. Lett. B439, 366-372 (1998). [hep-lat/9806026]. 

[17] http : //physics . indiana. edu/~sg/milc .html 

[18] A. Hasenfratz, F. Knechtli, Phys. Rev. D64, 034504 (2001). [hep-lat/0 103029]. 

[19] A. Bazavov etal. [MILC Collaboration], Phys. Rev. D 81, 114501 (2010) [arXiv: 1003.5695 
[hep-lat]]. 



-10- 



